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ORBIT DETERMINATION AND DIFFERENTIAL-DRAG CONTROL 
OF PLANET LABS CUBESAT CONSTELLATIONS 

Cyrus Foster* Henry HallamJ and James Mason^ 


We present methodology and mission results from orbit determination of Planet 
Labs nanosatellites and differential-drag control of their relative motion. Orbit 
determination (OD) is required on Planet Labs satellites to accurately predict the 
positioning of satellites during downlink passes and we present a scalable OD 
solution for large fleets of small satellites utilizing two-way ranging. In the second 
part of this paper, we present mission results from relative motion differential-drag 
control of a constellation of satellites deployed in the same orbit. 

INTRODUCTION 

Planet Labs Inc. is a purpose-driven space and data analytics company based in San Francisco, 
California. The company designs, builds and operates large fleets of Earth observation nanosatel¬ 
lites to image the entire planet at an unprecedented frequency. Planet Labs aims to provide universal 
access to information about the changing planet to enable both commercial and humanitarian appli¬ 
cations. 

Planet Labs’ concept of operations differs fundamentally from other optical remote sensing satel¬ 
lites in that the spacecraft always point their cameras towards nadir and are always on. Satellites 
always image while over land with reasonable solar illumination, and tasking is only done at the 
downlink stage. A single orbital plane of nadir-pointing satellites can image the whole Earth every 
single day if there are enough of them that their imaging swaths overlap. Lor the design altitude 
and telescope field of view of our instrument, a minimum of 105 evenly-spaced satellites is re¬ 
quired to ensure daily coverage as the Earth rotates under the orbital plane. Evenly spaced satellites 
additionally allow for the optimal use of the limited number of ground stations available for down¬ 
loading all of the imagery data. Eor these reasons, station keeping is strongly preferred over random 
distributions of satellites along the orbit. 

Operating fleets of satellites this size presents unique challenges, not least of which are orbit 
determination and station keeping control. The necessity to track each satellite accurately, and the 
ability to deploy, disperse, and then station-keep is critical to the optimal operation of the fleet. 
These capabilities were identified early on in the design of Planet Labs concept of operations, and 
the ability to station-keep without the complexity and expense of onboard propulsion drove us to 
develop in-house differential drag control technology. This paper presents a status update on the 
significant investment made over the last four years towards solving this problem. 
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On orbit determination, we present an introduetion to our UHF radio ranging eapability, and the 
resulting orbit determination eapability that allows us to traek our satellites independently of, and 
more aecurately than third party service such as JSpOCs. 

On differential drag, we present the basic physical concepts, the enabling attitude configurations, 
our constellation deployment strategy, and the methodology and mission results from differential- 
drag control of their relative motion. 

ORBIT DETERMINATION 

Orbit determination (OD) is required on Planet Labs satellites for multiple operational purposes 
including on-orbit attitude control, image capture and downlink scheduling, georectification of pay- 
load imagery, conjunction assessment, and most critically, pointing of ground station antennas dur¬ 
ing downlink passes. Imagery is downloaded from a satellite via a high-bandwidth X-band link 
requiring ground station antenna pointing errors of less than 0.2°. An ephemeris for the satellite 
must be known in advance of a downlink pass to compute an appropriate Two-line element set 
(TLE) or azimuth-elevation file for a specific pass. TLEs published by fhe Joinl Space Operafions 
Cenfer (JSpOC) are nol accurafe or reliable enough for nominal operafions, especially for safellifes 
in low orbifs (<450 km) where afmospheric and spacecrafl affifude uncerfainfies yield large perfur- 
bafions wifh orbif propagafion. JSpOC spacecrafl idenlificalion cross-lagging issues have also been 
encountered as a consequence of a large co-deployed heel of small spacecrafl. 

We presenl a scalable OD solulion for large heels of small safellifes. Two-way lime-of-highl 
ranging measuremenls are laken al UHE ground sfalions and are subsequenfly hi fo an orbif using a 
high-precision propagalor wifh a numerical force model. Orbifs are Ihen propagaled info fhe fulure 
lo produce a TEE or pass-specihc azimulh-elevalion hies. We presenl fhe melhodology for hfling 
an orbif lo Ihis dala and mission resulls. 

The safellifes, or Doves, fall info Iwo orbil-regimes: sun-synchronous inclinalion al 600 km alli- 
lude and sub-ISS orbifs (<400 km al 51.6°inchnalion). Al publicafion dale, Planel Eabs currenlly 
has 12 safellifes in fhe sun-synchronous orbif, and 25 safellifes in fhe sub-ISS orbif. 

Quality of JSpOC TLEs 

Especially when deployed in large numbers info low-allilude orbifs (e.g. via release from fhe 
Inlernalional Space Slalion), cubesafs can presenl a challenging problem^ fo JSpOC or any olher 
enlily aflempling lo hack Ihem via oplical and radar means. The combinalion of small radar cross- 
seclion and rapid evolulion in fhe orbif due lo varying afmospheric densily resulls in sparse, noisy 
measuremenls lhal musl be correlated and fagged by JSpOC’s semi-automated systems. Because 
each individual measurement provides no unambiguous identihcation of which object it represents, 
some cross-tagging inevitably occurs, leading to misidentihcation in the Satellite Catalog and/or 
corrupted solutions (when measurements from more than one object are mistakenly included in the 
same orbit ht). Accuracy of the TEEs published on Space Track* (as compared to truth from RE 
and laser ranging, GPS and phased-array imaging) has been highly variable and difficult to predict. 
Over the course of a few months, position accuracy of the TEEs issued for a particular object may 
vary from better than 1 km up to several hundred km of error. 


*https://www.space-track.org 
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Ranging data acquisition 

The aforementioned TLE quality issues, discovered early in Planet Labs’ Flock 1 mission, ac¬ 
celerated development of an independent means of orbit determination. Replacement firmware for 
the satellites’ UHL Low Speed Transceiver (LST) was developed and uploaded to the spacecraft 
already on-orbit, to allow two-way time-of-flight ranging measurements to be made. Unlike the 
high-bandwidth X-band link, the LST can tolerate ground station pointing errors in excess of 10°, 
allowing bootstrapping from an initially coarse orbit estimate. 

The LST ranging packets are tagged with the serial number of the satellite to be interrogated, en¬ 
suring positive identification and eliminating the possibility of cross-tagging. Ranging stations use 
commodity Yagi-Uda antennas and positioners, and can be deployed worldwide at very low cost. 
Since many satellites may be simultaneously in view from any site, a program was developed to 
automatically decide which satellite to interrogate based on various heuristics, including the satel¬ 
lites’ range and elevation, angular distances from current antenna position, distribution of recently 
gathered ranging data and operator-assigned priority factors. 



Figure 1. Example of ranging data acquired during a pass vs. JSpOC TLE 


Approximately 350,000 ranging measurements are gathered daily across three dedicated ranging 
stations. The Icr accuracy of each measurement is 650 m. Planet Labs recommends that all cubesat 
missions include some means of independent orbit determination, and to facilitate this has released 
an open-source version of the Low Speed Transceiver design. 

Orbit fitting to ranging data 

We fit an orbit to acquired ranging data using a high-precision propagator and least-squares solver. 
The propagator numerically integrates cartesian state vectors using the Adams-Moulton variable- 
order and variable-step method. We use a 10x10 EGM96 spherical harmonic gravity field and 
the NRL-MSISEOO atmosphere model with latest space weather and Earth Orientation Parameters 
(EOP) data from Celestrak*. 

To fit a state vector to ranging observations, we make iterative corrections to an initial guess 
using a batch Gauss-Newton least-squares solver as described in Algorithm 1.^ The initial guess 

*http://celestrak;.com 
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Algorithm 1 Orbit Determination Algorithm 
1: procedure Fit orbit given ranging observations 
2: get state veetor inital guess X from previous run 

3: eompute diagonal weight W from sensor 1-sigmas: Wi^i = — 

'^ranging 

4: while state veetor position update is large do 

5: propagate state veetor to epoehs of ranging observations 

6: eompute partials matrix A = and residuals Ay = Range^jigg^g^j—Rangepj.Qp 3 gjjjgj 

7: eompute eovarianee matrix P = {A"^WA)~^ 

8: adjust fit span and fit variables as appropriate based on P (on iteration only) 

9: update state veetor with 6X = PA^WAy 

10 : end while 

It: end procedure 


Table 1. State vector fit elements depend on quality of observations 

Nominal case Poor observability Very poor observability 
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is obtained from a state veetor produeed the last time it was solved for and propagated forward in 
time to the present epoeh, but an approximate guess from a JSpOC TLE ean also be used. While it 
should be possible to obtain an approximate state veetor from only the ranging data via Initial Orbit 
Determination (lOD), to date there has never been a need to do so. 

We deeide whieh elements of the state veetor to fit based on various eriteria for the eovarianee 
matrix P. For the large majority of eases, we ean fit the full state veetor X, but in some eases we 
have to resort to only fitting a subset of state veetor elements beeause of poor observability, sueh 
as when the satellite has reeently been deployed and may have only had a single ranging pass (see 
table 1). 

For example, if the magnitude of the position eovarianee is too high, we fit position and veloeity 
veetors eomponents in the orbit plane (R and S unit veetors in the RSW satellite frame^). If the 
position eovarianee is still too high, we then resort to fitting the satellite’s mean anomaly via an 
along-traek proxy. Similarly, the ballistie eoeffieient BC is fit if its eovarianee is small enough, its 
inelusion judged independently of the position and veloeity eomponents. We are also eareful not 
to use a fit span that is too long when the satellite is in a high-drag environment to ensure partial 
derivatives sensitivities do not exhibit highly non-linear properties. We fit the inverse BC beeause 
it is numerieally more eonvenient and exhibits more linear behaviour. 

Similarly, we ehoose an appropriate span of ranging data based on position and BC eovari- 
anees from P. Typieally, fit spans of 3-4 days are seleeted to be appropriate for the 600 km sun- 
synehronous orbit, while a shorter span of 0.5-1 day is ideal for the high-drag sub-ISS orbits. 




GPS data 


Some of the satellites earry an experimental software-defined GPS reeeiver, whieh due to power 
eonstraints and radio-frequeney interferenee eannot be operated eontinuously. However, the GPS 
ean provide sparse, highly preeise position/veloeity/time measurements against whieh ranging data 
and determined orbits ean be validated. 

To produee a high-quality ’’truth” ephemeris, we take segments of time that feature densely and 
eonsistently-aequired GPS observations. We then fit a state veetor (ineluding the ballistie eoeffi- 
eient) to all available GPS observations within an appropriate fit span, and sequentially repeat this 
process to cover the span of available GPS data. Finally, we evaluate position at desired times by 
interpolating between states propagated from bounding fit epochs. 

While the GPS data is not always available on the current generation of satellites, making it 
unsuitable for routine orbit determination, it is very useful as an independent source of observations. 
For example, the GPS data has allowed us to find and accounf for biases in fhe ranging dafa. We 
expecf GPS fo be fully-funcfional on fufure safellife builds, providing a more accurafe, if higher 
lafency, source of observafions for aufomafed operafional orbif determination. 

Orbit prediction results 

In figure 2, we show posifion offsefs for a ranging-derived ephemeris and JSpOC TLEs, as com¬ 
pared fo a high-qualify GPS-derived ephemeris. The ranging ephemeris is produced using all avail¬ 
able ranging dafa, and we inferpolafe befween bounding TLEs published during fhe period of inferesf 
fo produce a JSpOC-derived ephemeris. The comparison is made for safellife Elock lc-5 (NORAD 
CATID 40038) af reference epoch 2015-02-24 14:00:00 UTC, during a 12-day period where GPS 
measuremenfs were acquired regularly. Af fhe lime, fhe safellife was in a 590 x 620 km allilude 
sun-synchronous orbif. 

We observe a posifion offsel for ranging-derived and JSpOC TLEs of 0.45 and 1.45 km RMS 
respectively (maximums of 0.86 and 3.46 km) when compared fo fhe high-qualify GPS ephemeris 
in fhe 24 hours after fhe reference epoch. These resulls are represenlalive of whal we observe 
for Planel Labs salelliles in sun-synchronous orbils. Salelliles in fhe lower sub-ISS orbif will also 
fealure similar behaviour af limes bul can also exhibil higher posifion offsels in fhe 10-20 km range 
during periods of inconsislenl allilude maneuvers and space wealher variabilily lhal make orbif 
filling and propagalion wilh a fixed ballislic coefficienf less appropriate. These effecls are especially 
pronounced as fhe safellife encounters an increasingly dense almosphere as ils allilude decays. 

Eigure 3 shows posilions offsels for a ranging orbif solulion produced only using ranging dafa 
lhal was available before fhe reference epoch (fhe Iasi ranging poinl used in fhe til was acquired af 
2015-02-24 11:45:03 UTC). If is compared fo a high-qualify posl-processed GPS-derived ephemeris 
fo show fhe predictive capabilily of a lypical solulion from fhe orbif delerminalion rouline. Also 
shown is fhe posifion offsel befween fhe GPS-derived ephemeris and fhe lalesl JSpOC TEE lhal was 
available before fhe reference epoch, exemplifying fhe accuracy observed if using TLEs for orbif 
prediclion. 

Eor prediclions, we observe posifion offsels for ranging-derived and JSpOC TLEs of 0.56 and 
1.82 km RMS respectively (maximums of 1.03 and 3.84 km) when compared fo fhe posl-processed 
GPS-derived ephemeris. These posifion offsels are quoled for fhe 24 hour interval following fhe 
reference epoch, fhe period of grealesl inferesf for accuralely predicling X-band passes, and are 
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representative of on-orbit behaviour across the fleet. These results highlight the advantage of our 
in-house OD capability over reliance on JSpOC TLEs. 


5| ! ! - 

— Ranging post-process 
- - JSpOC post-process 
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Time since fit epoch (days) 


Figure 2. GPS Truth vs. Post-processed Ranging and JSpOC solutions 
(Post-processed uses all available data) 



Figure 3. GPS Truth vs. Predictions from Ranging and JSpOC 
(Predictions use only data available before epoch) 
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Responsible Operations 

Planet Labs’ mission operators, some having previously worked in this field at NASA, are acutely 
aware of the threat that debris poses to future satellite operations. We care deeply about it and are 
determined to operate openly and responsibly in space to ensure sustainable and universal access 
to low-Earth orbit, which we view as a limited and shared resource. We have put in place several 
best-practice policies to limit the negative impact of our satellites. 

Launch Low Selecting orbits that have low altitudes is the single biggest factor for reducing 
Planet’s impact on the space environment. Launching low reduces the time spent in orbit and ensures 
that we are flying below the congested 700-1000 km orbital band. In the unfortunate event of a 
collision, debris fragments will also have a short lifetime, further reducing the negative impact. 
Planet’s model of many small, relatively cheap, and short-lived satellites allows us to launch into 
orbits with lifetimes of a few years, as opposed to larger satellites that prefer to maximize time in 
orbit. We always ensure that our satellites are well below the UN 25 year guideline, and we will 
never launch a satellite into an orbit with a predicted lifetime greater than 25 years. 



Figure 4. Semi-major axis altitude of NORAD-tracked space objects in LEO. Planet 
Labs cubesats fly under the congested 700-1000 km region and exhibit lifetimes of less 
than 25 years. 


Open orbit data The greatest limitation with current collision predictions is the poor quality 
of TLE-derived orbits, which results in many false positives. To help increase the accuracy of 
collision predictions Planet maintains an orbit catalog of its own, which is produced using ranging 
data from several ground stations and is publicly published hourly to http: //ephemerides . 
planet-labs . com. We are a member of the Space Data Association (SDA), which streamlines 
the delivery of our orbit solutions to other operators. In the event of a predicted close approach, 
our high-precision ephemerides are directly available to the other operators to help them better 
understand the risk and decide on a course of action. Planet also reports all launches and orbits to 
JSpOC, and shares tracking data to assist them in identifying our satellites upon deployment. 
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Conjunction Response It is inevitable over time that high risk conjunctions will occur, and if 
the risk is deemed unacceptable we may want to take preventative actions. Planets Labs’ satel¬ 
lites are not equipped with propulsion, so avoidance options are limited. For all conjunctions with 
active, maneuverable satellites, we will contact the operator and point them to our public-facing 
ephemerides server to ensure that they have the best-available ephemeris for our satellite. If the con¬ 
junction is with any other object and has a small miss distance (and similarly small uncertainties), 
then we command the Planet satellites into a minimum interaction area attitude configuration at the 
time of closest approach. Since the collision probability is a function of the combined hardbody 
radii of the two conjuncting objects, this simple attitude maneuver can result in a 30% reduction in 
risk of collision with a 1 m-sized object, and the gains increase if the other object is smaller. It may 
also be feasible to use differential drag to avoid conjunctions, given enough lead time. We have 
developed a draft policy for differential-drag avoidance maneuvers, and are currently investigating 
the control authority and uncertainties involved before moving ahead with implementation. 

DIFFERENTIAL-DRAG CONTROL 

In the second part of this paper, we present methodology and mission results from relative motion 
differential-drag control of satellites deployed in the same orbit. 

Line scanner for the Earth 



Figure 5. Line scanner configuration for Planet Labs Earth imaging satellites 
Satellites are equally spaced in orbit and numbered to provide uninterrupted imaging 
coverage along the orbit path. The Earth rotates beneath, yielding global coverage on 
a daily basis. 


The ideal Planet Labs constellation is launched on a single rocket and placed into a 475 km sun 


synchronous orbit. Given perfect station-keeping, the number of satellites needed to achieve daily 
revisits is 105. If station keeping is not perfect, then there will be gaps in the coverage wherever 
satellites are spaced too far apart. This can be compensated for by additional satellites in the plane, 
to ensure some imaging overlap margin, although this needs to be accounted for in downlink bud¬ 
gets. In this way differential drag control authority can be traded directly against number of satellites 
(and amount of redundant data). To keep the cost and complexity of the fleet to a minimum it is 
strongly desired to minimize phasing control error. A constellation of evenly-spaced satellites has 
the additional benefit of minimizing ground station scheduling conflicts, and therefore maximizing 
the amount of imagery that can be downloaded from each spacecraft. 

Maneuvering without propulsion 

The Planet Labs satellites do not have propulsion systems, but since their orbits are low enough 
(<600 km) we are able to use differential-drag control to achieve the desired formation. Control 
is achieved by modulating the background attitude mode of the satellite when it is not imaging or 
communicating with a ground station. Since the satellites are highly non-spherical, different atti¬ 
tude modes yield different ballistic coefficients, different rates of atmospheric decay and therefore 
different rates of mean motion increase. By controlling the amount of time each satellite spends in 
a high-drag mode, one can ensure all satellites end up with the same mean motion resulting in zero 
relative speed. By adjusting when these high-drag maneuvers are made, one can target each satellite 
to a desired orbit slot relative to its neighbors. Unlike with propulsive thrusters, one is limited to 
controlling the along-track positioning of a satellite when depending on differential-drag control, 
since one is only effectively modulating the rate of decay due to atmospheric drag. 


a) 




Figure 6. Views of Planet Labs Dove satellite with cross-sectional areas 

a) Camera facing (200cto^) 

b) Solar-panel facing (1950cm^) 

c) Side-panel facing (370cm^) 


Due to operational considerations, we use two discrete attitude modes when a satellite is not 
occupied with imaging or downlink tasks: a nominal low-drag mode which is preferred for power 
collection, and a high-drag mode afforded by the satellite’s large deployed solar panels. The nominal 
low-drag and high-drag modes correspond to faces b) and c) respectively of figure 6 normal to the 
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velocity vector. These modes result in an approximate 5:1 ratio in surface areas exposed to the 
’’wind”. The control authority from this arrangement depends widely on altitude and atmospheric 
conditions, but ranges from ~ 1 km/day^ in a 600 km sun-synchronous orbit to ~ 50 km/day^ and 
greater in a sub-ISS 400 km orbit. 

Attitude maneuvers are achieved using the satellite’s Attitude Determination and Control System 
(ADCS), consisting of a star tracker and magnetometers for estimation, and reaction wheels and 
magnetic torquers for control. 

Control Methodology 

We define orbit slots at equally spaced angular distances from a fixed reference safellife labelled 
fhe ’’leader”, convenienfly chosen fo be fhe safellife in fhe lowesf orbif (highesf mean motion) since 
if is leasf likely fo require any high-drag maneuvers. We fhen create a fwo-componenf sfafe for each 
safellife as ifs mean Earfh-cenfered angular disfance fo ifs desired slof Oi, and ifs angular velocify 
relative fo fhe leader The objecfive of fhe confrol problem is fo fhen largel a final sfafe consisting 
of fhe desired orbifal slofs and zero relative angular velocity respecfively on all salellifes. 

Initial dispersion We assume fhaf fhe consfellafion salellifes are all launched fogefher and de¬ 
ployed wifh random AV componenfs relative fo fhe velocity vector. The laffer assumpfion ensures 
fhe initial separation of salellifes in fhe along-lrack direclion, as observed on fhe Plane! Labs Flock 
1-c launch from a Dnepr wifh ±lm/s along-lrack dispersion. If fhe initial relative velocities are nol 
sufficienf fo inilially disperse fhe salellifes in a timely manner, one can command salellifes fo each 
enler high-drag mode for an appropriale amounf of lime.^ 

Slot assignment The choice of which salellifes fo send fo which orbifal slofs relative fo fhe leader 
is made based on Iheir inilial dispersion. We can calculate fhe amounf of time required for each 
safellife fo reach each slof using componenfs of algorilhm 2, fhen choose fhe vector of unique slof 
assignmenls lhal minimizes lhal lolal amounf of lime to realize fhe formalion. The assignmenl of 
slofs needs to be performed after deploymenfs and inilial relalive velocities has been observed, and 
repealed if Ihere is significanl change to fhe availability of salellifes lhal perform differenlial drag. 

Algorithm We presenl in algorilhm 2 a procedure to issue high-drag windows fo a fleel of space- 
crafl fo achieve desired orbifal spacing. We define Ihese windows as scheduled periods of time when 
a spacecrafl has been commanded to mainlain a high-drag allilude when nol imaging fhe Earlh or 
engaged in a downlink pass. As a resulf of operafional limilalions, fhe confrol problem is non-linear 
since we only have Iwo discrele allilude modes to work wifh, a nominal low-drag mode and a ma¬ 
neuvering high-drag mode. The slrafegy is to firsl calculale fhe required amounf of high-drag time 
^ti,hd for each satellite fo malch ifs mean molion fo fhe leader (fhe safellife wifh fhe lowesf orbif 
which does nol need to perform high-drag). One fhen compules fhe amounf of time each safellife 
should wail in nominal low-drag mode Ati^wait before enlering Ihis high-drag window, such lhal 
each safellife malches fhe mean motion of fhe leader and enlers ifs desired orbifal slof simullane- 
ously. 

The predictions for fhe durafions Ati^hd, ^ti,wait are made wifh assumplions on fhe fulure slate 
of fhe almosphere and precise allilude of fhe spacecrafl. The accuracy of Ihis forecasting nalu- 
rally becomes worse wifh time, so we rerun fhe algorilhm on a regular cycle, and selecl an update 
frequency fo balance belween phasing accuracy, number of high-drag windows and availabilily of 
command upload opporlunilies. 

We presenl in figure 7 fhe simulaled behaviour of Iwo salellifes performing differential drag using 
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Algorithm 2 Differential Drag Control Algorithm 
1: procedure Create high-drag windows 
2: remove any previously issued high-drag windows 

3: compute 6i, 9i of each satellite i relative to previous leader (mean, not osculating) 

4: compute 0 of a satellite in high-drag vs. low-drag mode 

5: assign leader to be satellite with greatest 9 (i.e. lowest altitude) 

6: for each satellite i that is not the leader do 

7: offset 9i, 9i to be relative to new leader: 9i = 9i — 9ieader-: di = Ot — 9if,ader 

Or 

8: get required high-drag duration to null relative speed with leader: Ati^hd = — — 

^ 9 

9: get angle travelled during desired high-drag window: 9i^hd = -OAt'f^^ 

Or frd Or 

10: get wait time to target desired slot: Ati^rnait = —^^- 

Or 

11: create high-drag window for satellite i starting at tnow + Ati^wait for duration Ati^^d 

12: end for 

13: end procedure 


the procedure from algorithm 2. The two satellites initially remain in nominal low-drag mode due to 
their initial relative speed following deployment, then Sat2 (in a higher orbit) performs a high-drag 
maneuver to match mean motion with Sail (in the lower orbit), and timed to do so as the former 
enters the desired orbit slot at 100° from the latter. We distinguish the commissioning phase, from 
deployment to first attainment of desired slot, from the subsequent station-keeping phase, but the 
same control algorithm is used throughout. 



Figure 7. Simulated behaviour of two satellites performing differential drag. 

9 is the relative Earth-centered angle between the satellites. Thick lines represent 
high-drag windows 
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On-orbit results 


Flock 1-C orbital spacing relative to 090C vs. time 



Figure 8. Predicted Flock 1-C relative motion without differential drag. Note angles 
wrap around at 360°. Final state exhibits non-optimally placed satellites with various 
clumps and gaps, and high relative speeds. 



Figure 9. Actual Flock 1-C relative motion achieved by using differential drag, plot¬ 
ted from acquired ranging data. Thick Ones represent executed high-drag windows. 
Final state exhibits equaly spaced satellites with very low relative speeds. 
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We applied the differential drag eontrol algorithm to a lioek of 10 satellites in an approximately 
600 km altitude sun-synehronous orbit. Figure 8 shows a prediction of relative motion of these 
satellites if no differential drag control were implemented, while figure 9 shows the actual real- 
life behaviour of the flock from July 2014 through February 2015, as determined from ranging 
data. As one can see, we were successful in guiding the satellites to roughly equally-spaced slots 
and near-zero relative speed by using a sequence of high-drag windows over the course of six 
months. Satellites 0908 and 090C were a special case with assigned slots 2° apart, to demonstrate a 
relatively close formation, representative of the angular distance between equally-spaced satellites 
on upcoming launches. 

During this period, the flock was under the influence of algorithm 2, albeit with varying update 
frequencies due to higher priority operational constraints, resulting in periods without any high- 
drag windows across the fleet. Such delays implied that slot assignments had to be reassessed to 
ensure the time to realize the constellation was kept a minimum when differential-drag control was 
resumed. 

CONCLUSION 

From the very early stages of Planet Labs’ conception we have identified station keeping, and 
specifically ’’passive” differential-drag station keeping as a desirable capability. In line with Planet’s 
’’Agile Aerospace”, capabilities-driven engineering philosophy we have iteratively developed the 
in-house capability to perform both orbit determination and station keeping maneuvers. 

We have described and demonstrated successful results for orbit determination using two-way 
ranging data to reliably predict orbits with sufficient quality to perform X-band passes for data 
downlink from Planet Labs imaging satellites. While the current approach to orbit determination 
is sufficient for the high-altitude 600 km orbit, and most of the time in the low-altitude sub-ISS 
orbit, we recognize that further improvements can be made by incorporating short-scale attitude 
maneuvers in the orbit propagator, and operationalizing the use of GPS data. 

We have also detailed our approach to differential drag control and presented successful mission 
results from the phasing of 10 satellites to specific orbital slots and near-zero relative speed. The next 
challenge will be applying this approach to larger fleets of lOO-i- satellites deployed simultaneously 
in the same orbit, both for initial orbit phasing and to maintain equally spaced relative motion. 

Ranging and differential drag are capabilities central to the concept of operation of a large fleet 
of nadir-imaging satellites in a single orbital plane, designed to apply the line scanner strategy to 
capture imagery of the whole world every day. These techniques fit naturally into Planet Labs’ 
wider strategy for responsible and efficient operation of large fleets of propulsionless nanosatellites. 
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